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Abstract 

We develop a high-order energy method to prove asymptotic stabil- 
ity of flat steady surfaces for the Stefan problem with surface tension 
- also known as the Stefan problem with Gibbs-Thomson correction. 

1 Introduction 

The Stefan problem is one of the best known parabolic two-phase free bound- 
ary problems. It is a simple model of phase transitions in liquid-solid systems. 

Let Q C M" denote a domain that contains a liquid and a solid separated 
by an interface F. As the melting or cooling take place the boundary moves 
and we are naturally led to a free boundary problem. The unknowns are the 
temperatures of the liquid and the solid denoted respectively by f and v~ 
and the location of the interface F separating the two different phases. 

We shall assume that Q = T"~^ x [—1,1] where T^^^ stands for an {n — 1)- 
dimensional torus. Let us assume that the moving interface F(t) is a graph 
given by Xn = p{t,x'). Here p: [0,T] x f""^ — >R is some smooth function such 
that Uo<t<T-'^(^) T>0. Define the liquid/solid phase Q^(t) by set- 

ting = G x„^p(x',t)|. We note that fi = fi+(it)Ufi"(t). In 

order to formulate the problem we first specify the initial conditions. Let 
Fq = graph(po) be the initial position of the free boundary and : f2 — * R be 
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the initial temperature. The unknowns are the interface |r(t);t>o| and 

the temperature function v.Qx [0,T] ^R. We denote the normal velocity of 
r by V and normalize it to be positive if F is locally expanding Q~ (t) . The 
mean curvature of r(t) is given by 

.W=v.( , ). 
Vv/i + |Vp(()p' 

The Stefan problem with surface tension is now given by: 



dfV — Av = 


infi, t>0, 


(1.1) 


V = K{t) 


onr(t), t>0, 


(1.2) 


dnV = 


onT""ix{a;„ = ±l} 


(1.3) 


V=[d,v]t 


onr(t), t>0, 


(1.4) 


v{-,0)=vo 


inf2, 


(1.5) 


r(o) = ro. 




(1.6) 



Given v we write v"^ and v~ for the restriction of v to Q'^(t) and Q~{t), 
respectively. With this notation [(9j,f]l stands for the jump of the normal 
derivatives across the interface T{t), namely 

where z/ stands for the unit normal on the hypersurface r(t) with respect to 
Q^{t). If we replace the boundary condition fll.21) with 

v = on r(t), t>0, (1.7) 

then we are referring to the classical Stefan problem. 

The difficulties in dealing with the existence of solutions of the prob- 
lem (11.11) - (11. 6p arise from the nonlinear coupling between the temperature 
V and the boundary p. This connection is expressed through the boundary 
conditions (11.21) and (11.41) . The equation (II. 4p is a Neumann- type bound- 
ary condition for v. It is hyperbolic in nature as opposed to the parabolic 
diffusion process in the regions fl'^ and Q~ . 

From the technical point of view the first major obstacle for the analysis 
is the moving boundary. To deal with this issue we shall first transform the 
problem to the fixed domain by applying the so-called Hanzawa transform. 
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To this end let us fix a small positive constant « < | and choose a cut-off 
function 0gC°°(M) with Im((/)) C [0,1] and 

Define now a diffeomorphism 

and the function u{x' ,Xn,t) =v{Q{x' ,Xn,t)). Observe that u{x',0,t) = 
v{x' ,p{x' ,t),t) and at the outer boundaries dfl"^ ._ qpn-i x {x„ = ±1}, we have 
^|an± = u\gn±- This is the version of the transform first introduced by Han- 
zawa (cf. [I2])- In the new coordinates the heat operator dt — A transforms 
into a more complicated operator whose coefficients depend on the interface 
function p and the cut-off function 0. Following the calculations from j5j , we 
find that the Laplace operator A in the new coordinates takes the form 



AqU = A^fU + apUnn - ' V^j^'U^ - dpU 



where 



, 0Ap (0^)-|VpP , 0V(l + |0Vpn , . 

1 + 0'p (l + 0'p)2 (1 + 0'p)' ■ 
Furthermore, the operator 9^ in the new coordinates reads 

{dt)eu = dtU + epUn 

where 

0Pt 

ep:=- , , • 1-10 
1 + 0'p 

Note that RHS of (11.21) remains unchanged in the new coordinates. In order 
to transform the boundary condition (11. 4p into the new coordinates we first 
observe that 

{di)eu = diU- f^'^ dnU, l<z<n-l, idn)eu = , dnU. 

1 + p 1 + p 
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We thus conclude that at the boundary T" ^ x {x„ = 0}: 

d 

Vv\r = VeuU^=o = (V.m,0) - -^(0Vp, 1) = (V.m,0) -a„u(Vp, 1). 

The outward unit normal is given by uix't) = _ Tj-^^g ^]^q normal 

velocity V takes the form V = — ^=^==. Using the above expressions we derive 
the formula for ]j^|r. Namely on T"^-'^ x {x„ = 0} we have 

1 + 1 Vpp 



It is thus easy to see that the equation ( ll.4p transforms into 

dtp={l + \Vp\')[dnu]-. 
For the sake of notational simplicity we also set 

(p) := Vl + lVpP. 
The Stefan problem (11. ip - (11 .Gp now takes the following form: 

Ut - A^iU - apUnn + Bp ■ Vx'Un + CpUn = (1- H) 

u = n on T"-^ X {x„ = 0} (1.12) 

a„M = onT''-ix{x„ = ±l} (1.13) 

u{x,Q)=uq{x) xeQ (1-14) 

p{x' ,0) = po{x') x'er-^ (1.15) 

Pt = {pf K]; on T"-i X = 0}, (1.16) 

where we set Cp:=dp + ep with dp and Cp given by (II. 9p and (ll.lOp respectively. 
Recall that Op and Bp are given by (11.81) . In order to deal with the hyperbolic 
equation (I1.16P we introduce the regularization of the jump relation (I1.16p : 

Pt + eA^pt = {pf K]; on T"-i x {x„ = 0}. (1.17) 

We shall refer to the problem of finding the solution to (ll.lip - p.15p and 
(I1.17P as to the regularized Stefan problem. 

Notation. For notational simplicity we define for any multi-index p = 
(pi,...,p„_i) and sGN 

d^ = dtd^l...d>i::l. (1.18) 



Note that the operator acts only in directions tangential to the boundary 
qpn-i X {x„ = 0}. The Latin letters are always used to refer to the differen- 
tiation with respect to the time variable t and Greek letters to refer to the 
differentiation with respect to the first n — 1 spatial variables xi,. . . 
If each component of is not greater then that of /i and s' < s, we 
write {fi',s') <{fi,s). We write < if {fi' ,s') <{fi,s) and |yu'| < 

or s'<s. We also denote C^, =((^'5/))- For given functions u;:T"~^— *M 
and U:Q-^M., we denote u!i = dx^uj, z = l,...,n — 1 and Un = dxJJ- With 
X = (xi, . . . ,Xn) and x' = {xi,. . .,Xn-i) we set 

n— 1 n—1 
i,j=l i=l 

The Einstein summation convention is used throughout the paper when deal- 
ing with repeated indices. The letter C will stand for a generic constant that 
may change from line to line. 

We define the following high-order energy norms: 



|^j|+2s<2fe'^ ^ 

(1.19) 

\^i\+2s<2k•^ ^ 

|V^,a,^W(t)|2 + 2a^(t)|a,^V..W„(t)|'+(a^(t)9,^Wnn(t))'} 

+2 E / \dyuj,{t)\'{mr\ 



|/i|+2s<2fc 

where for given functions ujjifj: f"'"^ ^M, we define 



;i.2o) 



n-1 

,2/ ,\-3 



I^{V^u,V^tu) := I V'^l' i^y' - E (Vo^fc ■ VV^) (^)-' . (1.21) 



k=l 



Recall that a^(^t) is given by (11. 8p . It is crucial to observe that is a positive 
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definite bilinear form. Namely, 



n-l 



n IL J. 

Jjn-1 

(1.22) 

Note that we have used the Cauchy-Schwarz inequality in the last estimate. 
The instant energy S and the dissipation T> are respectively given by 

S = £{u,p):=£iu,p;p), (1.23) 

V = Viu,p):=Viu,p;p). (1.24) 

In the rest of the paper we shall always assume k>n, where n is the di- 
mension of the space the domain Q belongs to. Observe that the stationary 
solutions to the Stefan problem fll.lip - fll.161) are given by {u,p) = {0,p), 
where peM is a given constant. Note that S{u,p~ p) =S{u,p). The main 
result of the paper is the following theorem. 

Theorem 1.1 There exists a sufficiently small constant M* >0 such that if 
the initial data satisfy (mo,Po) satisfy 

S{uo,po)+ / Po- Mo(l + 0'po) 



<M*, 

then there exists a unique global solution to the Stefan problem U.ll\) - U.l()\) 
satisfying the global bound 

1 /■* 

£{u,p){t) + - V{u,p){T)dT<£{u,p){s), t>s>0. (1.25) 

2 J s 

Moreover, given the stationary solution {u,p) = {0,p) , such that 

P= Po- / Uo{l + (f)'po), 

then for such small initial datum there exist constants Ki,K2 >0 such that 
£{u,p){t) + \\p{t)-p\\l<Kie~^^^ for all t>0. 
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The proof of this theorem will strongly rely on careful examination of the 
regularized Stefan problem (11.111) - fll.lSp and fll.171) . For this purpose we in- 
troduce the appropriate energy norms incorporating the additional viscosity 
coefficient e. 

\^,\+2s<2k (1.26) 



+l4V^d^Auj,V^d^Auj) 
V,{U,uo-i)): = V{U,uo-i)) + 2t V / d"^ Auo\'^ . (1.27) 

\^l\+2s<2k'' ^ 

The above norms are the weighted versions of parabolic Sobolev norms given 
by 

ll(W,^)lk:= {ll^."W|lHHn) + PrV^llHi+e||5.^AV.;||^.} (1.28) 

|^i|+2s<2fc 

and 

||(W,^)|k:= {\\d':U,\\l,^^^ + \\dyU\\l,^^^^^^^ 

\^l\+2s<2k 

(1.29) 

Given ||'?/'||oo small enough so that (1 + 0^)^ ^ > ^-iid HV^/'Hoo bounded, 
we conclude that there exists C > so that 

^||(W,^)|k<^.<C||(W,a;)|k, ^m,u)\W<V,<C\\{U,uj)\\^^. 

In this sense the above norms are equivalent and this observation will be 
often implicitly used throughout the paper. The major part of the analysis 
will be concerned with proving the following result, which states that the 
regularized Stefan problem has unique global solutions with small initial data 
- independent of e. 

Theorem 1.2 There exists a sufficiently small constant M>0 independent 
of e, such that the following statement holds: if for given initial data (ug,pg) 
the inequality 



S.{ul,pl)+ / p^- / n^(l + M) 



<M 



holds, then there exists a unique global solution (m^,p^) to the regularized 
Stefan problem U.ll\) - U.15\) and ( [i.iTp . Moreover, 



1 f 

£,{u,p){t) + - V,{u,p){T)dr<S,{uo,Po), t>0. (1.30) 
^ Jo 

The Stefan problem has been studied in a variety of mathematical literature 
over the past century (see for instance ^1]). It has been known that classical 
Stefan problem admits unique global classical solutions in ([7j, [8] and 
|13j). Local classical solutions are established in [T2] and [IT]. 

If the diffusion equation (11. ip is replaced by the elliptic equation Av = 0, 
then the resulting problem is called the Hele-Shaw problem (or the quasi- 
stationary Stefan problem) with surface tension. Global solutions for the 
Hele-Shaw problem in two dimensions with small initial data have been es- 
tablished in [4J. In [2j, stability of the solutions close to the steady state 
sphere is established. Global stability for the one-phase Hele-Shaw problem 
is established in ^lOj. Local- in-time solutions in parabolic Holder spaces in 
arbitrary dimensions are established in [3]. 

As to the Stefan problem with surface tension, global weak existence 
theory (without uniqueness) is analyzed in [15] and [I9] . An existence theory 
is also developed in [1]. In [9] the authors consider the Stefan problem with 
small surface tension i.e. a-Cl if (II. 2p is substituted by v = aK(t). The local 
existence result for the Stefan problem is studied in [18]. In [5] the authors 
prove a local existence and uniqueness result in suitable Besov spaces, relying 
on the L^-regularity theory. 

We establish a global-in-time existence, uniqueness and exponential de- 
cay of classical solutions to the Stefan problem with surface tension near 
a flat steady state (Theorem II. ip . The major difficulty consists of proving 
Theorem 11.21 which establishes the existence and uniqueness result for the 
regularized Stefan problem with the energy estimate 




(1.31) 



where C does not depend on e. Combined with the smallness assumption on 
the instant energy £^ the estimate fll.3ip gives (11.300 and the global-in-time 
existence. For a fixed e we first construct local-in-time solution for the regu- 
larized Stefan problem (11.110 - (11.150 and (I1.17P . The crux of our method is 
the use of high order energy estimates, for the differential operator acts 
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only in tangential directions with respect to the boundary T"~^ x {x„ = 0}. 
This is very convenient when deriving the energy identities because the Neu- 
mann boundary operator commutes with d^. The diffusion equation (11.111) 
is then used to control high order derivatives of u with respect to the normal 
direction x„, as it is presented in Lemmas 13.51 and [3T6l We set up an iteration 
scheme, which generates a sequence of iterates {{u'^,p"^)}meN- Such iteration 
is well defined, but it breaks the natural energy setting due to lack of exact 
cancelations in the presence of the cross-terms. With fixed e, we crucially 
use the regularization to prove that p™)}meN is a Cauchy sequence in 
the energy space. As m — > oo the unpleasant cross-terms disappear and we 
recover fll.3ip in the limit. We conclude the proof of Theorem 11.11 by letting 

This work is the first step in our program of developing a robust energy 
method to investigate and characterize morphological stabilities/instabilities 
arising in numerous free boundary problems in applied PDE. In particular, 
in a forth-coming paper we are going to establish stability and instability(!) 
of steady spheres in the Stefan problem with surface tension. 

The article is organized as follows: In Chapter [2] we derive general energy 
identities for a model Stefan problem. In Chapter [3] the iteration scheme 
for proving the local existence is set up and the actual energy identities 
are derived, based on Chapter [2l Furthermore, some basic estimates are 
established, which are then used in Chapter H] to prove the crucial energy 
estimates. Chapter O is entirely devoted to the proof of the local-in-time 
existence nd uniqueness. The main results. Theorems 11.11 and II .21 are proved 
in Chapter [61 



2 Energy identities 

Let /=[0,g] for some 0<g<oo. The derivation of the energy identities 
crucially depends on the following model problem: 

Ut-A^>U-a^Unn = f on fix/ (2.32) 

U = Ax{^r^-^i^jXiA^r^ + G on T"-ix{x„ = 0}x J (2.33) 

dnU = on T"^^x{x„, = ±l}x/ (2.34) 

[Ur,]^={ut + eA'^tu){ijy^ + h on T"-^x{x„ = 0}x/ (2.35) 

We shall denote 
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For most of the identities we shall derive, only the leading term Ax (■?/')" 
in f l2.33p will be relevant. We can thus write the equation (12.331) in the 
alternative form 

W = Ax + g on T"-^ X {xn = 0} x /. (2.36) 

We define the energies and (for the model problem) by setting = in 
the definitions (I1.26P and fll.27p of and V^, respectively. 

Lemma 2.1 Let each of the functions U, uo, x and ip be five times continu- 
ously differentiahle with respect to the space variable and each of its spatial 
partial derivatives of order < 5 once continuously differentiable with respect 
to the time variable. The following identity holds: 

^aU,uj;tlj) + V,{U,uj;^)= [ {P + R}- [ {Q + S + T}, (2.37) 

at Jq Jjn-l 

where 

P = P{U,^, f):=fU~ {a^)nUnU, 

(2.38) 

R = R{U,ijJ):=f^-2f{B^-VMn)+Ul{a^)t-2UtUn{a^)t 
-2V^:Un ■ Vx'{a^)Un + 2A^./WW„(a^)„, 

(2.39) 

Q = Q{x,^,ij,g,h) = VuOf (Vx- V^) (^)-' + eAV^f (AVx- AVo;) (t/-)"^ 
-i{|V^P(^)7^ + e|AVa;|2(^)7^}+a;iVx-V((^)-^) 
+e/\VuJt ■ V((^)"') Ax - {ujt + t/\^ujt)g - {ipfhU\j^-. , 

(2.40) 

S = S{x,uj,^,g,h) ■.= 2VuJt ■ (Vxi - Vut) {^y' + 2eAVuJt ■ (AVxt - AVcu^) {^y' 
+2Vxt ■Vm-')uJt + 2eAWut ■ V( (^)-') Axt - 2 ^ + t/^^Ut) (Ax + gt) 
-2h{i)fUt\j^-i, 

(2.41) 

T = T{x,uj,^,G,h)=A{x,uj,ij,g,h)+B{x,uj,ij,g,h) + 2AG{uJt + eA^uJt) 
+2AW|^^n-l/^(V^)^ 

(2.42) 

Here, the functions A and B are given by: 

A := 2Acut (Ax - Acu) (V^)"' - 2AuJti'itlJj{xij-^^3) (V')"' - I V'cu|' 

+2uj,tVuJi ■ V( (^A)"') - 2VuJt ■ V{{'^py^)AuJ + uJjkOJ^k'^p3'^P^t {i^y'^ + (VcJfc ■ Vtp) ' {tpy 

-2u [i'ii'jUJkt {i'y^]i-i^ii^jUJikt I + 2ukt^ [i'ii'jUJij {ijy^]^-iJiijjUijk (V^) 

(2.43) 
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and 

B ■= 2eA^uJt (A\ - A^c^) {tpy^ - 2eA^uJttpiiJj {Axij - Auj,,) {tfj)' 
-e| V^Acup + 2eAcJitVAcJi ■ V((^)"^) 

-2eA Vwt ■ V((V')"') A^cj + eAu;JkAu;ikiJjiJ^t 0^) + e ( AVcUfc ■ V^) ' {tp); 
~2eAujjk^[iJjii'jAujkt {i'y'']i - i^ii^jAujikt {ipy 
+2tAuJkt[\iJii^jAuij - i)ii)jAuijk{ijy 

+2e{A(AxW-'-V^.V^,X.,W"')-(A\W"'-V^.^,Ax.,W"'^ 

(2.44) 

Proof. We start by multiplying the equation (12.321) by U and integrating 
over Vt. By a direct computation, 



(2.47) 



Jn Jn ^ Jt"-i Jn 

(2.45) 

where P{U,ip,f) is given by (I2.38p . Using the boundary conditions (12.361) 
and (I2.35p . we obtain 

- {i^f = -{u^tAx + eA'^.Ax {^y' } 

-{^fh[Ax{^r'+g]-iuj, + eA'uj,)g. ' 

Integrating by parts in the first term on the RHS of (I2.46P above we arrive 
at 

- [luAx {^r' + eA^uJtAx m-'} = 
f^„.,V{ut (^)-i) . Vx + eAVcu, ■ V( Ax (^)-') . 
By the product rule, the integrand on the RHS of ( I2.47P can be written as 

Vut{^)-'-Vx + eAVuJt-AVx{4^)'' + uJtVm~')-Vx + eAuJtVm-')-AVx. 

In each of the terms Voo {ip}^'^ -Vx and AVc^f AVx('^)~"'^ we set x = ^ + 
(x — uj) to obtain 

^9t||Va;|2(V;)"' + e|AVcu|2(V;)-'}+Vcur (Vx- Vcu) {tpy' + 

+eAVcur (AVx-AVcu) (^)-' - i| | Vcuj^ Ve| AVcup | + 
+UtVx ■ V((V^))-i + eAVa;i ■ V((^)-') Ax 

(2.48) 
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Thus plugging ([23HD into (1^:^71) yields 

- / {i^f[Un]-+U = \dt [ \\Vcu\m;)-' + e\AVcu\m;)-']+ [ Q{x,uj,^,g,h), 

(2.49) 

where Q{Xi^i''Pi9i^) is given by f l2.40p . To complete the derivation of f l2.37p . 
we take the square of the equation (12.321) and integrate over VL: 



-2 / {ijf[Un]-4At + 2 [ (^)'[W„];A,,W= / i?(W,^,/), 



(2.50) 



where R(U,'ip,f) is given by (12.391) . The goal is to evaluate the two integrals 
over T""^ on LHS of (12.501) using the boundary conditions (I2.36P and (I2.35p . 
We first treat the integral /j„_i (V')^ [^n]+Wt. Integrating by parts in the 
leading order term, we obtain 



-2 



(2.51) 



where SiXi^ii^iQih) is given by (I2.4ip . Note that the expression (12.410 is 
obtained similarly to (12.480 by setting x = + (x — cu) in the leading order 
terms. 

The second integral over T"~^ in the identity (I2.50p is more delicate. 
We shall make use of the boundary conditions (I2.33P and (I2.35P to evaluate 
it. The relation (12.330 is used is to exploit the full algebraic structure of 
the curvature-type term Ax{'ip)~'^ —ipiipjXij i"^)'^ ^ which is important in the 
energy estimates later on. We have: 

{^)'[Un]+A,'U = A{Ax{^r'-iJiiJJX^A^)~')i^t + eA'u;t) , . 
+AG{uJt + eA^uJt) + {ijfhA^,U. ^ ' ' 
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Observe that 

>-l ; ; / / \-3 



^ -1 JT"-i 

Jjn-l Jjn-l JT"-1 

2/" VcJi-V((^)"^)Acj + 2 /" AuJt{Ax-Aiu){ijy^- 

(2.53) 

Here, just like in (12.481) we substituted x = + (% — cu) in the leading order 
terms. Note that we have repeatedly used integration by parts. Integrating 
by parts twice, we obtain 

-2 / uJkkt'ipiiljjUJij{ilj)~^ = 2 / uJkt'ipi'ipjUJijk{ipy^ + 

Jfn-l I. 

-2 / Uikt^pi^pjUjjkW'^ -'2 / ujjkl [(pitpjUJkt{4^y^].-ipi'ipjUJikt{'ipy 
+2 / ujkt\[ipiipjUJij{ipy^]i^-ipiipjUJijk{ipy 

(2.54) 

We now single out the t-derivative in the first term on RHS of (12.541) to 
obtain 

n-l 

-2uikti^itpjUJjk {ipy^ = (^^fc ■ ^^i^YW'^ } +uJjkUJiki^jiJit (^)"^ 

k=l 

+ {Vuk-V^Y{^)f. 

(2.55) 

We combine the identities (12.531) . (I2.54p and (I2.55P to conclude 



-1 

dt \V'u\'{^)-'-dJy2 (V^fc-V^)'(^)-n+ / A 

=dt [ /^(vVvM+/ A 

Jfn-l Jj"-i- 
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where and A are given by fll.2ip and fl2.43p respectively. 

In the e-dependent part on RHS of fl2.52p we set u* = Au and x* = ^X- 
We can write 

We may now apply the same computation as in fl2.53p to conclude 

2e / A(AxW~'-^.^,X^, W"')AV = 9i / el4V^Auj,V'Au)+ [ B, 

where B is given by (12.441) . We combine the above identities to write the 
final form of the second integral over T""^ in the identity (12.501) : 

2/ {^f[Un]+A,'U = dt{[ L^{V^iu,V^Lu)+el4V^Aiuy^Auj)] 

Jjn-l I. Jjn-1 ) 

Jjn-l 

(2.56) 

where T is given by (12.420 . By summing the identities (I2.45P and (12. 50 p . 
plugging (ICTj) in ([23SD and fl^3T]) and (1^3^ into fl^3UD and collecting 
terms, we conclude the proof of the lemma. □ 



3 Iteration scheme and the basic estimates 

We shall set up an iterative scheme in order to solve the regularized Stefan 
problem locally-in-time. For given p"^ and Cauchy data UqEC°°{Q), pg G 
C°°{T^~^), we solve the following problem: 

<+i - A,,M™+i - Op^M^+i + Bp,n ■ V,'<+^ + Cp™M^+^ = (3.57) 
^m+i^^m on T"-^x{x„ = 0} (3.58) 

9„M™+i = on T"-ix{x„ = ±l} (3.59) 



Here 



u 



"^+\x,0)=ul{x), p^+\x',0)=pl{x'). (3.60) 



'^™-V-(^). (3.61) 
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The solution to the problem 03.571) - (13.601) exists and is smooth (see Chapter 
4 of [T^). Having obtained u"^~^^, we solve the equation 

p™+^+eAVr^^ = (p'")^[<+^]; onT"-ix{x„ = 0} (3.62) 

for p™'"^^. We aim for proving the convergence of the sequence {u"^,p'^) to 
the solution of the regularized Stefan problem in the energy space. Applying 
the tangential differential operator (recall (11.181) ) to the equations (I3.57p . 
f l338|) and fl332|) . we obtain 



where 

(3.63) 

^?;:.= E C^">r'Ap"^e"i^'((p"^)-^)+a,^(Vp--V((p-)-^)), (3.64) 

<ImI+s 

-{d^ipTpTpi] {pn'') -pfpTd'sPT, {pn'']. 
K,s= E c^''di{pr'+e^'pr')dn:{{pn-\ (3.66) 

Im'I+s' 
<ImI+s 

For any / G N let us define 

£':=£,{u\p'-p'-^), V':=Vlu\p'-J-^), (3.67) 

where and T)^ are defined by (11.261) and (11.271) respectively. Setting lA = 
d^u^+\ u; = d>!p"^+\ X = d^p"', iP = p"', f = fZ, 9 = 9Z^ G = GZ and h = 
the identity (12.371) implies 



d 
Jt 



£"^+\t)+V"'+\t)= [ {P"' + R"'}- [ {Q"' + S"' + T"'}. (3.68 
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Here P"" = J2M+2s<2kPZ and R"", QT , 5"^ and T"^ are defined analogously, 
whereby 

P^^=P{d^,u^^\p^Jl,), Ql=Q{d^,p-',d^,p^^\p'-,gl,,hl,) (3.69) 
and 

i?™ :=i?(9>"+\p'",/,":J, Sl,:=S{d^,p"^,d^,p"^'^\p^,gl,,hl,), 

The inequality (11.221) implies that the instant energy is positive definite. 
In order to estimate P™, i?™, Q™', 5™ and we first need to establish some 
basic auxiliary estimates. 

Lemma 3.1 The following identity holds 

dt{J^u'^-'\l + <P'pn}=dt{l^^ r^'}. (3.71) 

Proof. We multiply the equation (13.571) with (1 + 0'p™) and integrate over 
Q. We thus obtain 

/ (l + 0'p™)«+i-A.,n™+i)- / ^^^r^<r^' + ^ [ 0Vp™-V.mr^ 



n 

n 



n 



Integrating by parts we have J^—(f)p]^u'^^^ = J^cft'p^u"^^^. Using this 
identity, we obtain 

(l + 0'p™X+i+ / {l + <P'pne,u^+' = dJ [ (l + 0V™)n™+i 
Jn ^ Jn 

Observe that the integration by parts implies 

-(l + 0'p'")Ax'W™+'= / ^'V.'P'^-V.m'^+i 

Jn 

and /^0Ap"^<+i = -/^0Vp'^-V^/<+^ Thus 

- [ (l + 0'p™)A,m"^+i+ [ 20Vp"^-V,m™+^+ / 0Ap"<+i = O 
Jn Jn Jn 
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Note further that 

■l + |0Vp""p\ (02)'|Vp™|2 0"p'"(l+|0Vp""P) 



l + ^'p™ / 1 + 0'p™ (l + 0'p"')2 

Using integration by parts again we have 



1+0'p™ 7f^V l + 0'p" (l+0'p'")2 / " 

This finishes the proof of the lemma. □ 

The importance of this identity is reflected in the fact that it allows to 
control terms with purely temporal derivatives of p"*"*"^: 

Lemma 3.2 There exist positive constants K and 9q<1 such that for any 
9 <9q such that if 

I pl- I ul{l + c^'pl)\<9, 

k 

S'^<9, ||Vp"-^||oo<l and ^||c'pp'"||2 <irv^+^^, 

p=0 

then ||Vp"^||oo<l and EJ=o 1 1 V"^^ I b < KV£"'+^ + 9. 

Proof. Observe flrst that the assumption on p™" -1 implies that (p™"^) < ^2. 
Using the Sobolev inequality, we obtain 

1 1 VP*" 1 1 oo < a 1 1 VP*" 1 1 ^,(„+i)/2 < a V^S"" <aV29. 

Thus, choosing 9o< guarantees ||Vp™||oo<l- By (I3.7ip we have 

^T"-i Jo. Jn p=o\ P J 

Also, /^„_,p™+i = /^(u™+i + 0'M™+V'") + /Tr„-iP^- + Thus, for 
0<s< A; — 1, using the Cauchy-Schwarz inequality, definition fl3.67p of 



17 



and the main assumption in the statement of the lemma, we obtain 

s+l 



dsPt 



m+l 



< \\ds+iu"'^% + Cj2\\9pu"'^'\\LHn)\\ds+i-pP'' 
p=i 

k 



p=0 



Furthermore, 



m+l 



<C|hi™+l||2 + C|K+l||i.(0)||p'"||2 + 



Pi- / «^(l + 0'po) 



< C + C ^ I |<9pP"| I2 + 0. 

p=0 



By the Poincare inequahty we get 

k k 



J2\\9pP"'^%<cY.\\9p^P^^'\\^ + cY.\\ I (3-72) 

p=0 p=0 p=0 ^ 



The first term on the right-hand side is estimated by Cy£^^^, by the def- 
inition (13.671) of £^^^. By the previous two inequahties and the assump- 
tions of the lemma, we can estimate the second sum on RHS of (13.721) 
by C + C (kVS^ +e)+e. Keeping in mind that ^™<^, we 
choose 1<K large enough and 9 < small enough so that 



J2 1 \9pp""^^ 1 12 < kVs^+ e. 

p=0 

In the following we shall work under the standing assumption 

k 
p=0 



□ 



(3.73) 



with 9<9q where 9q is given as in Lemma 13.21 
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Lemma 3.3 Let ^ gC°°(J,M) and J CM. an interval such that every deriva- 
tive of ^ is uniformly hounded on J. 

(a) Let {fi,r) be a pair of indices such that \fi\ + 2r <2k and {\fi\,r) ^ (0,0) . 
Then there exists a positive constant C such that 

M[^{\^p""?)\\\m<C^ (3.74) 

and 

M[mP'^?)\\W^<^^- (3.75) 

(b) There exists a positive constant C such that 

\\d'^Vp^\\H^<CVv^, (3.76) 
where + 2s<2k. Furthermore, 

\\d^s+M'^P""?)\\2<CVv^ (3.77) 
for all {fi,s) satisfying \fi\ + 2s <2k. 

(c) For any pair of indices (/i,s) such that +2s <2k there exists a pos- 
itive constant C and a small parameter A such that 

1 1 d^pT'^^ 1 1 2 + 2e 1 1 Ap^+^ 1 1 ^ + 1 1 A Vr^^ 1 1 2 < y ^""^^ + C\V"'+^ . 

(3.78) 

Proof. Part (a): Let a = /i + r for any given multi-index of length n — 1 
satisfying |r|<l. Let first r = 0. By assumption |a|>l. Using Moser's 
inequality (cf. [11], Lemma 5) and Leibniz' rule (cf. [11], Lemma 4), we have 

ir[e(|Vp"^P)]||2<C max (||e(^)(|Vp"^p))|UI|Vp'"||H-i)p°(|Vp™p)||2 

<CP°(|Vp'"P)||2<C||a"Vp"^||2||Vp™||oo<CV£^. 
Let r>l. 

r 

5.ie(ivp"r)]=9"{^ c,^^'\\'^pn')dsA\^pn')--dsA\'^pn')} 

(1=1 siH HS(i='" 

r 

= E E E ^<^C7,,..,,^,9J^H|Vp-r)...5J^^(|Vp™r)9^-^[e^'^)(|Vp-r)]. 

d=l siH 71 H 

Sj>0 —a 

19 



For any z = l,...,d, we have d]ii\Vp^\') = j:tLoEs<,Ci,sdfVp^d]iZ;Vp"'. 
Thus if |7j| + 2si<k + l, then by the Sobolev inequahty 

and analogously | Vp"^| U < Cv^, implying \\d2i{\V p'^\'^)\\oo<C8'^ . 
If |7d+i| + we use the Sobolev and Moser's inequality to conclude 

||^7d+l [^W(|Vp"^n]||oo<^^V^. 

If there exists l<j<c? such that |7j| +2sj > fc + 1, then |7j|+2sj<fc + l for 
1 <z < (i, i 7^ j and additionally, |7d+i|</c + l. Thus we can estimate the term 
containing 7j in superscript in L^-norm and the remaining terms in L°°-norm. 
If on the other hand |7i | + 2s j <k + l for every 1 < i < ci, we estimate the term 
[^'^'^^(l Vp™"!^)] in L^-norm and the remaining terms in L°°-norm. We 
conclude that 

for the specified range of a-s and r-s. The inequality (13. 751) is proved similarly. 

Part (h): By ( l338ll . Ap"^ = m™+i (p"^) +prpf p^^ (p"')"^ Let 7 = p + r 
where |r| = 1. Applying 5J to the above identity, we get 

7',s' 

E c:i^:si:r,;^:dipTd2i^^^^^ 

— 7+s 

Observe that = since = = V ■ (Vp*" (p*")"^) on T""^ Let 

us fix (7',s') < (7,s). If |7'|>0 note that /^(9J, u"^^ = 0. We use the trace 
inequality and then the Poincare inequality on VL to deduce 

If |7'| = by the Poincare inequality and the trace inequality: 

||9.'«"+'||l2(t-i) < C||a,,V,,M'"+i|U2(Trn-i) 
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By part (a), \\d]Z], ((p"'))||2 < C + Cv^- Furthermore, if for some l</<4 
we have \'yi\ + 2si>k, we estimate the term containing 7; in superscript in 
L^-norm and the remaining terms in L°°-norm. Using part (a) we deduce 

\KpTdT,PTdT,p7]d]t{{n~')\\2<c{sn'^' m^pnin^- 

\fj,\+2s<2k 

Thus, summing over all pairs {p,s) yields 

J2 \\d]Ap"'\\2<{C + CVS^')Vv^ + C{S"'f/^ ll'^s^^P"'!^^- 

(3.79) 



\^l,\+2s<2k \fi\+2s<2k 

7=M+i",|t|=1 



Since Jj„-i Vp™ = and I7I > 1 elliptic regularity implies 

Y p^'vp"^iih^<c \\9:^pn\2. (3.80) 



\lJ.\+2s<2k \^\+2s<2k 

7 — M+T, I T I —1 



Combining fl3.79p . (13.801) and choosing sufficiently small we deduce the 
claim. 

Part (c): We apply the differential operator 9^ to the 'jump relation' (13.621) 
and take squares on both sides to obtain 

\\d>:pr'\\i+2ew,Apr'\\i+e'\\d>:/\'pr^^^^ 

Next 

Ws ( {pn' [«r 1 12 = I E^'r^r' ((p'")')e"i^' [uT%\\2 

fi',s' 

<c Y m'iipn')\U\d'::^'[ur%\\2+ 

m'|+2s' 

<fc 

+c Y (ipn")\\2m:^' 



m+1] 



+ 1 |L°°(T"-1) 



m'I+2s' 

>fc 



A A 

Here we assume S"^ < 1. Terms involving L°°-norm are first estimated by the 
Sobolev inequality. Then we use the standard trace inequality | |f | |l2(9q) < 
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xl I'^^l |L2(f7) + A| I Vf I |L2(f^) to bound the terms involving Observe that 

the same proof is easily adapted to yield the bound 



dspT^^\\2<CVv^, for s<k. 



(3.81) 



□ 



Remark. The estimate ( ]3.78p will play a crucial role in the energy estimates 
for the problem of local existence. For any pair (/i,s), + 2s < 2/c and > 1, 
the elliptic regularity and the estimate (13.781) imply 



I s rt II - _)^^2 



(3.82) 



Lemma 3.4 There exists a positive constant C such that 

(a) Forr>0, \fi\ + 2s<2k and (|/i|,s) ^ (0,0) 

\\d^dnr{ap^)\\mn)<CVS^ 

(b) For r>0 and \p\ + 2s<2k 

||a,^5„.(i?p-.)|U2(^)<Cv^ 

(c) For r>0 and \fi\+2s <2k 



mdnrM\\mn)<cvs^+cVv^. 

Furthermore, for r>0, \fi\+2s <2k ~1 

md^r{cp^)\\L2(^n)<CVS^. 

Proof. Part (a): We note: 

d^dnria,^) = 9,^9„.((l + |0Vp'"n(l + 0V-)-2) 

r 

= d^(j2dn.F,{\<pVp"'\)d^r-.F2{<P'n 

p=0 



(3.83) 



(3.84) 



(3.85) 



(3.86) 
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where := 1 + and -F2(x) := (1 +x) ^. Observe that 

r—p 

d=l s\ + ---^d 

— r—p 

and thus 

r—p 

r—p 
d=l fi',s' 

where Gd{x) : = x'^. By Moser's inequahty, for any pair of indices {y,q) such 
that \u\+2q<2k, ^ (0,0) we have 1 [[F2^'^^((/)'p™)] | Il^cq) < C'V^ if 

|z/|>0, and ||(9^[[F2^'^^(0V'")] ||L2(f^) <Cv^+^ if |z/|=0. We have used 
Lemma [321 Similarly, if |z/|+2g<2A;, we have 1 1^^^ [^^(p™)] | ^2(^5) < 
when |z/|>0, and ||9g^[Grf(p™)] |U2(f^) <Cv^+^ when |i/|=0. Thus 

11^^'' [[i^f (0'p'")]e"i^' [Gd{pn]\Wm<Cy/£^, (3.87) 

where we hit the terms with lower order derivatives with L°°-norms, de- 
pending on whether |p'|+2s'<A; or \p'\ + 2s' >k. Additionally, we use the 
assumption that 9 <1 and S^<1. This implies \\d^dnr-p 
CVS"^ for 0<p<r and |p|+2s<2A;. In the same way we prove 
\\d^dnpFi{\(l)Vp"'\)\\L^^n)<CV£^ for 0<p<r, \p\ + 2s<2k and i\pls)^ 
(0,0). Using the same idea of estimating lower order terms with L°°-norms 
as in the proof of ( ]3.87p . we conclude the proof of part (a). The proof of part 
(b) follows in a completely analogous way. 

Part (c): To prove part (c), recall that Cp^n = dpm + epm, where dpm and Cpm 
are given by 

^ 0Ap™ _(02)'|vp™p 0V'"(i+i(/)Vp'"n ^ 0p™ 

i+0'p™~ (i+0'p'")2 ^ (i+0'p"^)3 ' 

The analogous proof as in the part (a) imphes that | |9^c)„r (dpm) | |^2(q-) < 
C\/£"^. Furthermore, since \\d^pT\\<CVv^ if |p| + 2s<2fc, we use the 



(0'p-)(p^ 



m\d 
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same method as in part (a) to prove \\d^dn'-{epm)\\i^2(^Q-^ <CVV^. We thus 
conclude fl3.85p . From the definition of S"^ (cf. f l3.67p ) and the assumption 
of Lemma [321 we have | |a,Vrl I < for +2s<2fc-l, |/i|>0; also 

I |<9sP™| I <C^/£^ + 6 for all s<k — l. Now we use the same method as in the 
proof of part (a) to deduce (13.861) . □ 



Lemma 3.5 For all pairs of indices (/i,s) such that \fi\+2s + r <k + ^ + 3, 
the inequality | |9,^<9„rM™+i| |i2(f^) < Cv^l^ holds. 

Proof. We prove the claim by induction in r. In case r = l the claim is 
obvious from the definition of 8"^. Let the claim be true for all r<d for 
some 1 < -i? < fc + n/2 + 3. We have to prove the claim for r = 'd + l, i.e. 

I I |z.2(f,) < C^/8^\ (3.88) 

if + 2s < A; + 1 — ^9 + 2. Let (/i,s) be such a pair of indices. Then 

=5^5Zc(/i',s',w)|9,^>„.Mr+i-9,^;9„.A,.m™+i+ 

Ijb' ,s' u)=0 

+ Y.c{^l\s\p)^^J,^nv{Bp^)^^^,z^:: dn^-,+.v,m'^+^+ 

pii^s" p=0 

(3.89) 

Observe that \\d^', dn^.uT^^\\L'^i^^)<C^/£^ and 1 1(9,^,'5„™ A^/m"^+1| |L2(n) < 
(^y/^m+i fQj, g^j-^y triple of indices (/i',s',w) < (/ijS,-;? — 1) by in- 
ductive hypothesis. Further, by Lemma 13.41 and the Sobolev in- 
equality I 1 9^,, dnpBpm I |L°°{n) <C\/8'^, Wd'^ii dnvCpm\\i^ca(yi-^ < C\fE^ and 

||afri;X^-i-»(a;i)||L-(n)<Cv^+C, for (/i",s",p)<(/i',s',w)< 
(/XjSji? — 1). Applying these estimates to the above identity and using 
the inductive assumption, we obtain 
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where we recall the smallness assumption on S"^, specially S^<1. This 
finishes the proof of the lemma. □ 



Lemma 3.6 // + r + 2s < 2A; + 1, 

mdnrU"'^' I |i2(f,) < CVS^. (3.90) 

Proof. We prove the claim by induction in r. In case r = l the claim is 
obvious from the definition of S"^. Let the claim be true for all r<?9 for 
some 1 <2k + l. We have to prove the claim for r = i) + l, i.e. 

if \fi\ + 2s + <2k. Let (/i,s) be such a pair of indices. Then 

w 

fj,"^s" p=0 
w 

(3.91) 

We analyze separately the case when \^'\+w + 2s' <k and \fi'\+w + 2s' >k. 
Case 1. In the case \^'\+w + 2s' >k note that 

\fi~IJ,'\+2{s-s') + ^-w<2k-k = k. 



By the Sobolev inequality and Lemma 13. 4[ we have 
If (/i",s",p)<(/i',s',w), |/i"|+P + 2s"<A; and 



25 



where we have used the Sobolev inequahty and Lemma [23] respectively. Sim- 
ilarly 1 19^„ dnP {Bprn^ I iLoo(f^) < C\fE^. This implies 



I d^; ) 1 2 1 {c^,:^: v.m"^+^ ) 1 2 {d^z^; {a-i ) ) 

where we have used the inductive hypothesis to deduce 
Analogously 

2 



If on the other hand +p + 2s" > /c, we use the Sobolev inequality and 
Lemma [3.51 to get 

where we note that 

|/i'-/i"| + (w;-p+l)+n/2 + l + 2(s'-s")<A; + n/2 + l, 
so that Lemma 13.51 is applicable. In analogous fashion it follows 



We also note that 



l^ 1 12 

L°°{n) 



Observe that we have used the inductive hypothesis in the last inequality 
above. This completes the first case. 
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Case 2. In the case \^'\+w + 2s' <k, by the Sobolev inequahty and 
Lemmas 13.41 and 13.51 

and \\d^J,dnv {cp^)d^',:^',' dnr.-P+iu"'^^ Wloo^^) < CVS^ for (/i", s",p) < 
{li\s\w). By Lemma [331 Ws-s' dn^-^-^{%^)\\L^{n)<CV£^ + C. We also 
note 

By the Sobolev inequality and Lemma 13. 5^ 

We combine the above estimates to conclude | |(9^9„.j+im™'+^| |i2(Q) 

and this completes the second case and finishes the proof of the lemma. □ 



4 Energy estimates 



Lemma 4.1 Let K and 9<9q he given as in Lemma \3. 2[ There exists 0< 
L<6 and such that if 



JT"-i Jn 

and for some m G N 

sup £"'{t)+ [ V"'{T)dT<L 
0<t<T^ Jo 



L 
< — 
- 2 



then 



p=0 

sup S"'+\t)+ / V"'+\T)dT<L. 
0<t<T'- Jo 
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Proof. With the preparation from Chapter [21 we are ready to estimate RHS 
of fl3.68p term by term. Note that the assumptions of Lemma 1X2] are fulfilled 
and we are thus able to use Lemmas 13.31 - 13.61 in the forthcoming estimates. 
Let be an arbitrary, but fixed pair of indices satisfying \^\ + 2s<2k. 

Term J^PJ^^s- Recall that Pl[^, = P{d^u'^+\p"'J"') is given by fl339|) . 
where P is given by fl2.38p and by fl3.63p . Thus, combining (12.381) 
and fl3.63p we can estimate the first term on RHS of /^^-P^ 



n 

<c 



m+l^ 
nn I 



m+1 



m+1 



c y: +c 

\^l'\+2s'>k 



\lj.'\+2s'<k 
(|m'1,s')7^(0,0) 



(4.92) 

In the first sum, observe that \\d'^, apm\\j^cx,(^Q) <C^/S^ for \jj,'\ + 2s' <k, by 
the Sobolev inequality and Lemma 13.31 In the second sum observe that 
I Idglg/ u^n^ I < CvP™+^, by the Sobolev inequality and Lemma [3161 By 
glancing at (12.380 and (13.630 the second term in the expression PJ^^ is given 
by f^d^{Bpm ■'Vx'U^^^)d^u"^^^. It is estimated in a completely analogous 
way and is bounded by Cy/S^'D'^'^^ again. By (13.690 . the third term on 
RHS (I2.38P renders the third term in JqPJPs- We have 



d^{cp^u^+')d^u^+'\<CY,\ f d^:cp^d^:^'u^+'d^u^+' 



=^ E +^ E 

\fj,'\+2s'<k \fj.'\+2s'>k 



(4.93) 



In the first sum we estimate \ \d^Cpm \\^ like above ( since \fi'\+2s' <k). In the 
second sum, for \fi'\+2s'>k, \\d'^I^' u';^+^\\Loo^n) <CVt^ by the Sobolev 
inequality and Lemma [331 By Lemma [373| \ \d'^, Cpm \ \i^2 <C\jD^. Finally, the 
fourth term of f^PJ^g (again use (I2.38P to identify the fourth term and the 
equations (I3.69P and (13.631) to plug in the appropriate values), is estimated 
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by using the Cauchy-Schwarz inequahty 



fl 


<l 







m+l I 



(4.94) 

Term /j„-i<5™s" '^^^ proceed with the estimates for the expression 
/Tn-i<5".' where Q™, is given by fICTD . where Q is defined by (E^QD, c/;;;, is 
given by fl3.64p and h^,, is given by (13.661) . The first two terms of jf^-^^li^s 
are the cross-terms and they deserve special attention. For any r7>0, 

<r]\\ d^Vp"'^^ 1 1 2 + - f 1 1 <9,^ Vp"* 11^+11 dt'Vp""-^^ \\i)< r/P'"+^ + -(£"" + 8'^+^). 

(4.95) 

e / dyAp]^+' ■ I A9,^Vp™ (p'")"' - A9,^Vp"+' (p'")"' ) 

Jjn-1 I J 

<r/6||a,^AVpr'||^ + — (||a,^AVp"||^+P,^AVp^'"+^||^) (4.96) 



C 



Observe that the constant C does not depend on e. The third _term in 
Ijn-iQ'^^s is given by (CTD and (IXB^]) . Note that 9t((p™)"^) 



1\ _ 2Vp"'Vp^ 



By Lemma [331 we conclude || (p'")^ ^ ||oo <CV"^. We then obtain 



'jpri — 1 



< CD'" (I |9,^Vpr+^ 1 12 + e| |a,^AVp'"+^ 1 1^) < 



(4.97) 



To estimate the fourth term of /]p„_i Q^^ (which is obtained as the fourth term 
of (12.401) together with the definition (13.691) ). we use the Cauchy-Schwarz 
inequahty to get 



/ a,vr^9,^vp'"-v((p"^)-^) <||W+^||2P,^Vp^'"||2||v((p"^)"^) 



< 



(4.98) 
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Analogously, the fifth term in Jj„-i Q^^ is estimated as follows: 



/ 



AVa,Vr'-V((p'")-^)Aa> 



< 



\\^feAV^^pT+^\\2\\V{{p'^y^)\\oo\\^feA^^p'^\\2<CVV"^+^Vv^V^ 



■)m+l^ 



(4.99) 

We first note that the sixth term of (12.401) contains g. Note that in fj„-iQ^^s 
g = g^^si where g"^,, is defined by (13.641) . We shall first estimate || (7^5 112 and 
||y^V(y'^g||2 and then use the Cauchy-Schwarz inequality. For any |/i'|+s'< 
|/i| +s, we have 

The inequality follows by estimating the term with smaller order space- 
derivatives in L°°-norm, which can then be estimated by the Sobolev in- 
equality and Lemma IX^ Similarly, recalling (I3.67p : 

\\^ev{di/\p^d^y^:{{p^)-'))\\2+\\^fev{d':{vp^-v{{p^^^ 

Therefore Wg'^^^Wi+W^ftW g'^^,\\2<C\/S^\/V^ and we can bound the sixth 
term in j^^-iQ'^^s by 



JT"-i 



+ 1 1 A V pr"^^ 1 1 2 1 1 V g';^J\2<CVV'^+^ VS^ VV^ < C Vs^ (V"" + 1)'"+^ ) . 

(4.100) 

The last term in jj„-iQ^s is extracted from the last term of (12.401) . which 
contains h. By (13.691) . h = h^^ where /i™^ is given by (13.661) . For the nota- 
tional simplicity, we set /i^^ =: /ii + e/i2, where 

(4.101) 



<|m|+s <|mI+s 



Note that for \p' \ + s' < \p\ +s, we have 

\\{pn'd^:pr'd^:,^'{{pn-')\\2 ...r.^. 

<ll(p"')'llooPr'Pr''9,^ri;'((pn"')ll2<c'v^v/^^. ^ ■ ^ 



30 



Here, if \^i'\ + 2s' <k tlien \\d'^', pT^^\\oo<CV£^ and if \^i'\ + 2s'>k tfien 
ll^^ri^ ((p™)"^)||oo <C\J E"^. We use tlie Sobolev inequality and Lemma [3T3] 
to conclude the estimate. This implies \ \h^\\2<C^f£^^^pV™^ . Note that in 
similar fashion, for |yu'|+s'< + 



In other words | |\/e/i2| I2 <C\/ £''^\/T>''^^^ . From the proof of part (b) of 
LemmaEJl we deduce | I2 < CVV^ and also 1 1 v^9,^k™| I2 < CVV^ 

(recall here (I3.6ip ). Thus the last term of ^^n-iQ^^s bounded by 



{P^'K^sd's^"^ <C^||/^l||2||9,^«:'"||2 + C||v^/.2||2||v^9^^ 

'jpn — 1 



|2 



m+l 



(4.103) 



Term J^K^.,: Note that i?™^ is given by (13.701) where Z^^, is defined 
by (13.631) . Our first task is to estimate the first term of f^K^^, namely 
UfZ)'- Observe that 

7^(0,0) 

(4.104) 

If \fL'\+2s'<k, by Lemma E31 

I |<9^/ Opm I Iloo(s^) + I |(9^, Bpn. I |lcx,(Q) + | \d^, Cpm \ |ioo(Q) < C . 

Thus for |/i'|+2s'<A; and (/i',s') 7^ (0,0), RHS of flCTij) is bounded 
by If \fx'\+2s'>k, then + 2(s - s') < 1 and from 

Lemma 13.61 

In addition to this, for such (/i',s'), we use Lemma [33] to conclude 
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Thus, for every A > 



Analogously, 



for all < (yU,s) satisfying \fi'\ + 2s' > k. Combining the estimates for 

lyu'l + 2s' <k and + 2s' > k, we obtain 



Jn 



mcm+l 



(4.105) 



The second term of the expression R^^ (the second term on RHS of (12.391) 
and (I3.70p ) is bounded by: 

/ (a>r')'K-)t <ll(ap™)tllL^(n)P>r'||i.(^)<CV^P-+^ (4.106) 
Jo. 



"p"' Jn 



<CV^™I?™+i. Note 



and similarly the third term 

further that for the fourth term in J^R^^s (use f l2.39p and (I3.70p ). by 
Lemma 13.31 



(4.107) 

By Lemma [331 the last term in (last term on RHS of (IXBOj) and (KW ) 

is bounded by: 



<||A.,9>'"+i|U2(f,)||(a,™)„||i.o(f,)||9>r'llL^(f^)<C'v^^^™+'- 



(4.108) 



Term Jjn~iS^y. Note that 5*™^ is given by fl3.70p where S is given 
by flCTjl . c/;;;, by 'dMlD and /i™, by (13:661) . The first two terms on RHS 
of (12.411) are the cross-terms and in order to estimate them we shall exploit 
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the part (c) of Lemma 13. 3[ It turns out that the constants on the right-hand 
side will depend on e. Note that for any r], A > 



-1 



'J'n — l 

C,,^ m+l||2 I ^\\f)fJ'\7^rn\\2^ ^ _^ ^J^pm+l 



V 



(4.109) 



In the last estimate we have used the estimate (I3.82p . Similarly, 



c 



-V 



rn+l 



(4.110) 

The third term of J^^-i S'^s is given by the third term on RHS of 02.411) 
and flX7D]l : 



d^p^'d^vpT-viipn 



^71 — 1 



< 



m+1 1 



ll^^p; 



m+l^ 



Similarly, the fourth term in j^^-i SJ^^ (given by the fourth term on RHS 
of dMH) and (Km ) is bounded by: 



/ Adi^pTAvd^pr'-viipn-': 



< 



\\V~eAd^pT\\2\\V~eAVd^p, 



m+l I 



I V ( (p" ) " ^ ) 1 1 oo < C v^Vp^v^ 

<cV^(p™+r'™+^). 

(4.112) 

Note that the fifth term in Jjn-iS"^^,, by (ICTD and fl3T0|) ) involves the 
function (yf^^,, where g"^,. is given by fl3.64p . The crucial step in estimating 

is 



this term in JJ„-^S'^, is to observe that | |(^;^Jt| I2 < CV^VP^. This is 
proved by first differentiating g^,, with respect to t, and then in each product 
estimating the terms with lower order space derivatives in L°°-norm and the 
other one in L^-norm. The same method applies to show | |^/eV(5'^5)f 1 12 < 



CVE^VV^. Also observe that 

\\Adi:p^{pn;'\\2<\m 



p. m I 



{p"');^\\oo<cVs^Vv^. 
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Analogous proof shows that 

1 1 v^V ( A9,V™ (p™)7' ) 1 12 < Cv^VP^. 
Using the above inequahties and the Cauchy-Schwarz inequahty we estabhsh 

(4.113) 

Integrating by parts and using the analogous argument as in (14.1131) we get 



t \\2 



Jjn-l \ 

< 1 1 ^ed^Avpr' 1 1 2 1 1 v^v [Ad^p"" {pn;'+i 

(4.114) 

The sixth and the last term in /rj.„_i S"^ ( given through the last term on RHS 
of (12.411) and (I3.70p ) involves the function where h^,, is given by (I3.66p . 
Since d^^^!^ = V-d^ {V p"" {p"");^ ) , we can integrate by parts to obtain 



{pn'hZd^f^T- 



'^{{p"'fhz)-{vp"^{pn~'), 



We split h'^^^ = hi + eh2 as in fICTTD . For \p'\+ s' <\p\+ s 

\\'^{{pn"h)\\2<\M{pn')\\oom:pr^^^^ 
+1 1 (p-^)' I loo (p^; vpr ^e^i;' ((p-)~^) I b + p^^/pr ^e".^^ 

(4.115) 

Here we have used L'^ — L°° estimated by separating the cases \p'\ + 2s' <k 
and \a'\ + 2s' > k. Since 



||v^V((p")'9,^;AVriCi^'((p'")"'))||2<C'v^v/^^, 



we deduce 



\\V^Vi{p"'f h2)\\2<CV£^VV"^+\ (4.116) 

where /12 is given by flCTTD . Similarly, | |(9,\i (Vp*" (p'")"^) 1 12 < Cv^v^ 
and also 1 1 Ve9,\i (Vp™ (p™)"^) 1 12 < Cv^v^. 
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In summary, 



+ ||v/i((p-)^/.2)||2||v^5r(Vp™(p™)-^)||2 



(4.117) 



Term J^^^iT^^: Recall that T™^ is defined by f l3.70p where T is given 
by C;^, by' fl335|) and /i™, by (^M^. In particular the term A- the 

first term on RHS of (12.421) is given by (I2.43p . The first two terms of the 
expression A are the cross-terms. Using integration by parts, we obtain 



[ 9,^ Apr+' {d^^p"^ - d^Ap"^-^') {p^^y 

Jjn-l 

I 9,^Vpr+'-v((9,^Ap'"-9,^Ap™+^) {p"")-' 



m+li 1 2 
t I I2 



c 



(4.118) 



+ -||V(9,^Ap-(p-)-^)||^+ 



<A||9,^Vp 

^1 1 V(a,^Ap"^+i (p™)"') 1 12 < AI?'"+' + -^(^"^ + £™+^), 
A eA 

where we note that by the definition of S"^, | |<9^(Vp'" (p*")^^) 
for |p| + 2s<2fc. Similarly, for the second cross term: 



<XD 



m+l 



eA 



m I c'm+l^ 



(4.119) 



The proof of (14.1191) relies on the same idea as above; we first integrate by 
parts and then establish the estimate 



1 1 V [pTpj {d^sPT, -d>:pT'){pn-')\\l<^ (^"^ + ^ 

By A — (crossterms) we denote the sum of all the remaining terms in the 
expression A (recall (12.431) and the fact that ip = p"^ in our case) . Terms of 
the form (p™)"^, (p™')^'^, p™' and p^ for 1 < i,j < n — 1 are bounded in 

L°°-norm by Cy/S™-, by Lemma [3l3l Terms of the form (p™)"'^ are estimated 
by 1 in L°°-norm. Note that in the last two terms in the expression A (12.431) 
the leading order derivatives cancel out after the the product rule has been 
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applied within the parentheses. Using these observations and applying the 
Cauchy-Schwarz inequality, we conclude 



I A - (crossterms) | < C\f£^ (P™ + V 



m+l^ 



(4.120) 



Recall now that the term B (the second expression on RHS of (12.421) ) is given 
by (12.441) . The first two terms in the expression B are again the cross-terms. 
By part (c) of Lemma 13. 3^ we obtain 



]pn — 1 



< 



e^l |9,^AVr^' 1 12 + C| |9,^A V"! I2 + C| l^s'^^V"'^' 1 12 < 



(4.121) 



< ^^"+1 + (A + C^™)D™+i 
A 



C 



Analogously, we establish 



"Tpn — 1 



>-3 



< 



m+1 



{X + C8"')V 



c 



(4.122) 



We denote the sum of the remaining terms in the expression B hy B — 
(crossterms). The same idea as in the estimates for \ A— (crossterms) | works. 
It is important to note that we have canceling of the highest order derivatives 
within the parentheses in the last three expressions on RHS of (12.441) . In 
addition to that, we factorize e = ^/l x By the Cauchy-Schwarz inequality, 



I B - (crossterms) | < C\f£^ (P™ + V 



m+1 ^ 



(4.123) 



The third term of j^^-iTJ^g is given by the third term on RHS of (12.420 
together with (I3.70p . Recall that is given by (13.651) . In order to estimate 
it, we first integrate by parts. 



VGZs ■ d^Vp^' - 6 / A VG":, ■ d^AVp^ 



m+1 



The crucial observation is 1 1 VG™J I2 < Cv^v^, | |v^AVG™,| I2 < 
C \J E^\jT>'^ . Both inequalities follow in the standard way, by using — l? 
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type estimates and fl3.67l) . The third term of fj^-iTJ^^ is then bounded by: 

+ 1 1 v^A VG™, 1 12 1 1 V^d^AVpT^^ 1 12 < cVs^Vv^Vv^< cVs^iv"" + . 

(4.124) 

We integrate by parts in the fourth and the last term of /j-n-i TJ^^ (given by 
the last term on RHS of (KWf and (l3T0D ). Recall AU\jn-i = Ad{u 



m+l I 



Adi^HJ^. 



-[ Vi{p'^fhi)-dyK"'-e[ V{{p"'fh2)-dyK"'. 
By (I4.10ip . we set = /ii +e/i2- By the trace inequality, 



By M.llSp . we obtain 

[ Vi{p"'fh)-d^\/K"' <||V((p")'/li)||2||9,^V/t'"||2<Cv^I)"+\ 

Jjn-1 

(4.125) 

On the other hand, by Lemma 13.31 and L°° — type estimates, 
\\^d^VK'^\\2<CVV^. By ( Kim . 

e^_^V((p'")^/i2)-9,^V«:'"|<||v^V((p™)^/i2)||2||v^9,^V«:"^||2 ^^_^26) 



where we recall ( 14.1011) again. From the estimates (14.1251) and (I4.126p . the 
last term in Jj„-iTJ^^ is bounded by 

/ {fi^fKA^s^^^ <CV£^{V"' + V"'+^). (4.127) 

Jjn-l 

Using the identity (I3.68P and summing the estimates (I4.92p - (14.1140 
and (I4.117P - (14.1270 to get a bound on the right-hand side of the iden- 
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tity fl3.68p . we arrive at 

+C(ri + A)r''"+i + + £"'+^) + ™+i + ryp'" 

r? A 

T^Ae^ Ae 

(4.128) 

Now integrate in time over the interval [0,t] to get 

Jo ^ S<t VJg ^ Jq J 

s<t Jo Jo 

11 Cf Cf 

+Ct{- + -) supf™+i(s) + — sup^'"(s) + -— sup£™+i(s) + 



CXe 



V Jo 



[ V'^+\T)dT + 7] [ V"'{T)dT+^fsnpS'^{s) + snpS 

Jo Jo ^ s<t s<t 



+ SUp£™+^(s) [ V"'{T)dT. 
s<t Jo 

Note that by Cauchy-Schwarz inequahty we have 
sup 

s<t Jo 



(4.129) 



<sup yj£^+^{s)sj V"'{s)ds^ V"^+^{s)ds (4.130) 
<{[ r''"(r)rfr)^^^(sup^™+^(s)+ /" r'™+i(r)dr). 

Jo s<t Jo 

By assumption supg<^ f^"^ (s) + Jg'D'"(r)(ir < L, and thus from (14.1291) 
and Kim , for any t'~<t: 

Jo 2 7] Xe 

(sup^-+^(s)+ fv^^\T)dT)\L'/' + C{v + X)+Cti- + h + -^^ + ^ + ^ 
s<t Jo r] X r]Xe^ r] Ae 
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Since the above inequality holds for any t' < t, we obtain 

(snpS"^+\s) + fv"^+' (r) rfr) 1 1 - (l'/' + C(r/ + A) + + ^) + + ^ + ^ + 
s<t Jo 'IV 7] X r^Ae" r] Ae / 

<^ + CL'/' + -L + r^L+^. 
2 T] Xe 

We first choose r] small and then A small so that C{X + r]) + ^^ is small. 
Further, we choose t {t depends on e) and L small so that 

LV2 + C(r, + A) + Ct(l + l) + ^ + ^ + ^ + L<l 
1] X TjXe^ Tj Xe 3 

and 

^ + C(Lp + -L + r/L+^<^L, L<e. 
2 Tj Ae 3 

With such a choice of L and t=:T^, we obtain sup5<^e£^"*+^(s) + 
'D'^^^{T)dT < L and this finishes the proof of Lemma [4. 1[ □ 



5 Regularized Stefan problem. 

The principal goal of this section is the following local existence theorem: 

Theorem 5.1 For any sufficiently small L>0 there exists t'^>0 depending 
on L and e such that if for given initial data (mq,Pq) 



£,{ulpl-pl)+ / pl- / M^(l + 0'p^) 



ipn — 1 



L 
< — 
- 2 



then there exists a unique solution {u'^,p'^) to the regularized Stefan prob- 
lem U.ll\) - U.15\) and ( [i.!?] ) defined on the time interval [0,^^]. Moreover, 

sup £,{u%p';p')it)+ [ V,{u\p'-p'){T)dT<L 
o<t<t^ Jo 

and Se{u'',p''){-) is continuous on [0,t^[. 

Remark. Note that the constant L is independent of e. 

Proof. Convergence. Combining Lemmas 13.21 and 14.11 we obtain a 
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uniform-in-m bound on the sequence \^{u"^,p"^)^^. Our goal is to show that 
{(■u™,p™)}^ is a Cauchy sequence in the energy space. For any / gN let 
v'''^^ : = u'''^^ — u'' and 0"'+^ = p'+^ — p'. By subtracting two consecutive equa- 
tions in the iteration process, we obtain 

(5.131) 

V =A(r {p } 

= Aa"'{p"')-'-pTpfa'^{p'")-' + Gl, on r'-ix{x„ = 0}, ^ ' ^ 
5„t;™+i = on T"-^x{xn = ±l}, 

(5.133) 

=(^r'+eAVr^)(p™)-V/i:;, on T-'x{x^ = 0}. 

(5.134) 

Here 

/;:, = -5,.-V.,.;r'-V.t;r' + <„(a,™-a,™-0 

VTn / r> H) \ '^^ { \ yO . LOO j 

Gi, = Ap-i ( (p™)-^ - (p™-i> ) + p™-i {aTpf {pn-'+pr'^T (p"^)"' + 
prvr'((p™)"'-(p™"'»)' 

^"l ^"T- / „'TJ \ — 1 I 

9m — '~PiPj^ij\P I +tJm5 

= - (I Vp-^p - 1 Vp'-T) K]; (p"^)-' . (5.136) 

After applying the differential operator 9^ to the equations (15.1311) . (15.1321) 
and (15.1341) and singling out the leading-order terms we arrive at: 

d>:vr'-A^,v^^'-a,^v^^' = (5.137) 

(5.138) 
(5.139) 

where 

/; = 5."/^+ (5s"K'"W„^')-«P'"Wn''), (5.140) 
9'm = di:g°m+ E ^''Aa™e"i^'((p™)-^), (5.141) 

<ImI+s 
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ir,';r;: (5.142) 
{d>: [pTpTaT, ipn'')- pTpid's<, {pn~' ) , 

hL = d^h:^+ E d'::{crr'+eA'crr')d^:s^'{{p"'r')- (5.143) 



Im'I + s' 
<ImI+s 



As a next step, we use the identities from Chapter [2] to obtain the energy 
identities for the problem (I5.137P - (I5.139p . Respecting the notations of 
Chapter [2] we set for any I G N 



I 

(5.144) 



f = fl,g = g[,G = G'i,h = h[,U = d^v'^\ u = x = and ^ = p 

Additionally, we introduce the notations 

e': = £,{v',a'-p'-'), ■. = V,{v\a'-p'-'), 

where and are defined by (11.261) and (11.271) respectively. Using (12.371) 
and (15.1441) . we arrive at 

dt Jo, Jt"-i 

where 

pm,= J2 P{d^v"^^\p^J'J, r-:= E Rid>"''-\P"'JL) 

|/i|+2s<2fc \fi\+2s<2k 

q^:= E" Q{d>:a^.d^s<y'^^\p"',ChL). 

|^i|+2s<2fc 

5-:= E" S{d^,<T^,d^,<T^^\p^,g'M. 

|/^|+2s<2fc 

r:= E nd>'^,d':a'-^\p^,g'M- 

|^i|+2s<2fc 

(5.146) 

Here P, Q, i?, ^ and Tare defined by (EJHI), (ICTD . (E39D, (ESD and (EJZ]) 
respectively. Our aim is to prove that for suitably small t<t{e) there exists 
a A < 1 such that 



e"'+\t)+ I rf™+i(r)rfr<A(e'"(t)+ / d™(r)rfr). 
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We shall accomplish this by estimating the terms p*", r"^, g™, s*" and t™' on 
RHS of the identity fl5.145p . These estimates will be largely analogous to the 
estimates from Chapter |H However, due to the formally new terms d^fl^, 
d^g°^, d^G°^, d^h°^ appearing in the definitions fISTHOD . flSTHTj) . fl5ll2D 
and fl5.143p of g'^, G'^ and h'^ respectively, we need to make several 
preparatory steps. First, for any / gN 

I \0'p^ — 0'p^-A\hi{q.) + I \Bpm — -Bpm-i < C\ |cr™||//i+i(f^), (5.147) 

and 

I |cp™ — Cpm-i ||j^!(f^) < C||(J™||/^; + C| |cr™||j7i+2(f2). (5.148) 
Note that supo<«(E ^™(t) + jl V'^{T)dT<L. In particular, for |/i| + 2s<2A; 

||5>"^||/fMf^)'P>p-IU^(f^)'IIW"IU^(f^)<^v^- (5-149) 

Furthermore, \\d^Cpm\\^2(^^-^ <C\/L for \^\+2s <2k — l. We can now use the 
Sobolev inequality to bound the lower order derivatives of u"^, apm, Bpm and 
Cpm in L°°-norm by Cy/I. The major step is to provide the analogues of part 
(c) of Lemma [3.31 for the function cr"^+^ instead of and Lemma [3.51 for 
the function t)™+i instead of u™"*"^. By the boundary condition fl5.139p and 
the proof of Lemma [3731 part (c), we deduce 

I l^^r"^^ 1 1 2 + 2e 1 1 Act™+i 1 1 2 + 1 1 9,^ A 1 1 ^ < Ac/'"+^ + y ^ 

(5.150) 

and 

P^r^'l I2 + 2e| \d^Aar^'\\l + e^] |9,^AVr+' \\l< C + e'^V"' . (5.151) 
As in Lemma 13. 5[ 

I Id^dnrv""^' I |i2(f^) < Ce'"+^ + Ce'". (5. 152) 

The proof of 05.1520 is completely analogous to the proof of Lemma 13.51 
whereby, due to the addition of the formally new term d^f^ in the definition 
of we need to exploit the relations fl5.147p and fl5.148p . which are respon- 
sible for the occurrence of the term Ce™ on RHS of (15.1520 . We proceed fully 
analogously to the energy estimates in Chapter HI to estimate the right-hand 
side of the energy identity (15.1450 . The terms involving d^f^ and d^h'^ re- 
quire an additional care. In the estimate for the term J^p"^ = J^fm^s'^^^^^y 
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where is given by fl5.14UI) . we single out the term J^d^f^d^v'^'^^. Here 
is given by fl5.135p . Writing 

we estimate the lower-order terms in both sums in L°°-norm and the higher- 
order terms in L^-norm. By the Cauchy-Schwarz inequality and (15.1521) . 

I [ a,^(cpmw;r+^)a,^w'^+^|<C'v^y^^(v^+y^^). (5.154) 
Jn 

Similarly, | J^d^iBp^^ ■ V^'vl!^+^)d^v"'+^\ < Cy/£^{e"'+^ + (1^^+^). In analogous 
fashion, we find 



Jn 

(5.155) 

where we rely on (I5.147P and (I5.148p . The term r*" is of the form r"^ = 
(/m)^ + (rest), (r"^ is defined in (15.1460 ). The formally new term to estimate 
has the form J^{d^fm)^, where is given by (I5.135p . By the same splitting 
idea as in (I5.153P and the estimate (15. 1520 . 



Recall now that and t"^ are defined by (15.1460 . The last term in each 
of the expressions fICTD and fl2:i2|) has the form {p'^f h'^v^^'^^ljn-i and 
{p"^)'^ h'j^Ax'v"^~^'^\jn-i, respectively. By (I5.143p . the formally new term (with 
respect to h"^,, defined by (I3.66P ) is d^h"^. By (15.1360 . we conclude 

5.^/^;;.=-Ec'"'^"\[<];)e"^((ivp"^p-ivp-T)(p™)"')= E + E 

p',s' ii'+2s'<k fi'+2s'>k 

Hitting the lower order terms with L°°-norm and the higher order terms with 
L^-norm and using the trace inequality to estimate [m™] + ||2, we obtain 

C 



\\d^h°J\2<{XVV^+-VS^)V^, \\d^h°J\2<CV^^V^. (5.156) 

A 



43 



Note that v]^^'^\jn.~i = (Kpm — Kpm-i)t. It is easy to check that 

||a>r^^|Trn-l||2<C(||VVr||2+||V<||2 + ||VV™||2 + ||V(T"^||2). 

By the definitions of dr+^ and e™+\ \\di^v'^^\n-i\\2<^^V^+C^. 
Combining this with (15.1561) and the Cauchy-Schwarz inequahty, we obtain 

Choosing A = y/e we arrive at 



TTn — 1 

(J Q 



(5.157) 



We want to estimate the time integral of the right-hand side of the inequal- 
ity (15.1571) . For any t<t'',A>0, we have 

/ Vv^\f^VdrdT<\ ! d"'{T)dT + ^ sup e'^is) [ V"'{T)dT. (5.158) 
Jo Jo A 0<s<t Jo 

Similarly, 

/ Vv^e"'dT< f 6""+ sup e™(s) / V"^{r)dT<t sup e™(s)+L sup e"^(s). 

Jo Jo 0<s<t Jo 0<s<t 0<s<t 

(5.159) 

Since S"^ < L, for any A > 0, we obtain 

[ -VS^V?^Vd^dT<X [ d™(r)rfr + ^Lt sup e"^(s), (5.160) 
Jo ^ Jo Ae o<s<t 

and similarly, 

[ -^V£^e"'dT<-^tVL sup e"\s). (5.161) 
Jo Ve o<s<t 



Noting that ||A^/f 



m+l I 



equality and the estimate (15.1561) to get 
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2< -^ve™-, we can use the Cauchy-Schwarz in- 

C 



^n—l 



^/e 



(5.162) 



We get 

[ -^v^e'"dr<-t sup e'^(s) + L sup e"(s). (5.163) 

Jo ^ 0<s<t 0<s<t 

In analogy to the estimate fl4.128p together with fl5.154p . fl5.155l) . 05.1571) 
and ( I5.162p . we obtain 

d^+i^^m+i < (j^m^l ^ C ^ ^ C'e"^+i(i + 1 + ^ + ^) + 
at r) Xe r] X rjXe.* Xe 

C r— ^ C ^ C 



Integrating the above inequality over [0,t] and proceeding as in (14.1291) . us- 
ing fl57[58|l . fl5J[59|) . flSTTBOj) . fl516T]) . flSTTMD . we conclude 

e-+i(t)+ /" c;-+i(r)c/r<Ct(- + -^) sup e'"(s)+r7 /" rf'"(r)t;r+ 
Jo V Ae o<s<t Jo 

C^(- + T + 4^ + T^) sup e"^+i(s)+C(v^ + r/ + A + ^) / d"'+\r)dT 
T] X TjXe'^ Xe o<s<t V Jq 

+CX{ [ d"'{T)dT+ [ d'^-^\T)dT) + 



{^ + L + t+^t+^^t+-t) sup e"'is) + Cit + L + VL) sup e"^+^(s) 
A Xe y/e e o<s<t o<s<t 

(5.164) 

We choose r], A, L and 1=:^ < small, so that ri + CX<x<l 

^ A CL ^ ^ CL ^ cVl C ^ 1 

^^^;^+A^) + - + ^ + ^ + A?^ +7^ <^<5' 

Ct^(i + i + ^4)+C(f + L + v/L)<,<i 
and C{\/L + ri + X + ^^) <x< |- Taking supremum over [0,^"^], we arrive at 

sup e'^+i(r)+ / d"'+\s)ds<x\ sup e'^(r)+ / rf"^(s)rfs| + 



0<r<t^ 



+x\ sup e'"+^(r)+ /" d"'-^\s)]ds. 



45 



Therefore, 

sup e™+^(r)+ / (r'-\s)ds<-^{ sup e'"(r)+ / (r{s)ds\. 



0<T<t^ 

(5.165) 

We observe now that the conservation law fl3.7ip and the fact that a"^~^^{0) = 
v"'+^{0) = imply 

a™+^= [ t;™+i(l + 0'p"')+ / 
T"-i Jn Jn 

By the Poincare inequality, previous identity and the uniform bounds on p"^ 
and u"^ we get 



|a™+i||2<C||V(T'"+i||^ + C|| / a 



"n—l 



12 



m+l||2 1 /^\\ m+l\\2 i /^\\„,m\\2 ll^™||2 (5.166) 



<C||V(T™+^||^ + C||t;'"+^||^ + C||u'"||i2(^)||a'"||^ 
< ^6"+^ + la"*! 1^ < Ce"+^ + Id'"! 1^. 

With L and x so small that CL + 2C-r^ =: A < 1, we obtain by 05.1651) 
sup ||cr'"+'(r)||2<C^^| sup e'"(r)+ / ^"^(5)^4 + 
+CL sup ||cT"(r)||2<-| sup e'"(r)+ /" rf™(s)rfs| + ^ sup ||cr'"(r 



1 12 
2- 



Adding { supo<T-<te 6™+^ (r) + (i™+^(s)(is} to the both sides of the above 
inequality and using fl5.165p again we get 

sup |e™+^(r)+ /" c/"+i(s)rfs + ||a"+^(r)||2|< 

0<T<t^ ^ Jo ' 

C^^i sup e™(r)+ / (r{s)ds\ + -{ sup e™(r)+ / d'^{s)ds\ + 

1 — X'-0<r<t'! J 2 lo<T<t^ Jo ' 

sup A||a™(r)||^<A sup |e™(r)+ /" rf'"(s)rfs+ | |a'"(r)| j^j. 

0<T<t^ 0<T<t« '-JO 

Define the Banach spaces 

"{v,a)\\T),<oo 



X:=\{v,a) |(t;,a)|k + ||a||^<oo|, Y:={{v,a) 
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where IK-jOlb^ defined by fll.28p and fll.29p respectively. 

Let Z ■.= L°° [X; [0,^]) ©L^ (Y; [0,t']) . Since A can be chosen arbitrarily small, 
we have proven that the sequence ((^',0"'))^^^^ satisfies ||(f"*'''^,cr'"^^)||z < 
A'lK^^ja"^)!!^, for some A'<1. This implies that ((M',pO);gfs} is a Cauchy 
sequence and converges strongly in Z. Thus the whole sequence ((M',pO)jgN 
converges to a solution {u^^p'^) of the regularized Stefan problem in the orig- 
inal energy space. In addition to this, passing to a limit in (13.711) . we obtain 
the conservation law 

dt[j^u{l + (t>'p)]=dt[j^ ^p}. (5.167) 

Uniqueness. We want to prove uniqueness in the class of functions 
(u,p) satisfying s\^iQ^^^^t£e{u,p){t) + V^{T)dT<L, where L may be cho- 
sen smaller if necessary. Let us assume that there exists another solution 
(f,(j) satisfying the same initial conditions {v{x,0),a{x' ,0)) = (uo{x),po{x')) 
and the bound supQ^f^^e£^{v,a)(t) + V^{v,a){T)dT < L. After subtracting 
them and setting w: = u — v, T:=p — a, we obtain 

Wt - Aa;/W - ttpWnn = f*, 

(5.168) 

w = Ar{p)-' + g* = Ar{p)-'-pip,Tij{p)-' + G* on T"-ix{x„ = 0}, 

(5.169) 

[wn]^={n + eA\){p)-^ + h* on T^~'x{xn = 0}, (5.170) 
where 

/* = -Bp ■ V^'Wn - CpWn + (ftp - a„)Vnn + [Bp - B^)V^'Vn + (c^ - Cp)Vn, 

G* = Aa{{py^ - {cxy^)+aij [npj {py^ + aiTj (p)" VcXiCTj ( (p)~^ - (a)"^)) , 

9* = -piPjnj{py^ + G* 

/^* = -(|Vpp-|Vap)K];(p)-^ 

We use Chapter [2] to derive the accompanying energy identities. To this end 
we setU = w, ip = P, = X = ''': f = f* : 9 = 9* y G = G* and h = h*. Here p takes 
the role of p™^^ and a the role of p"* and additionally, the cross-terms vanish 
since a; = x = r. With k>n sufficiently large, the regularity assumptions of 
Lemma 12.11 are fulfilled. We are thus naturally led to the following energy 
quantities: 

£* ■.= £,{w,t;p), V* ■. = V{w,t;p). 
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In addition to this we define P* = P{w,p, f*) and analogously Q*, R*, S* and 
T*. Using the identity fl2.37l) . we obtain 

±e-+V=f{p- + R-}-r {q'^S'^T-}. (5,171) 

Our goal at this stage is to prove the inequality of the form 

j^£* (t) + V* it) < C£* (t) + cVlV* (t) , (5. 172) 

which would enable us to absorb the multiple of T>* on the right-hand side 
into the left-hand side and then use the Gronwall's inequality to conclude 
that S*{t) = for any t>0. It is essential that the constant C in the above 
estimate does not depend on e so that the smallness bound on L remains 
independent of e. That the identity fl5.172p indeed holds, follows analogously 
to the energy estimates from the Chapter H] applied to the right-hand side 
of fl5.17ip . Here we strongly exploit the uniform bounds on £^{u,p) and 
£e(f,cr). In particular we know that 

'j''p)t||oo; 1 1 ('3'p) 1 1 oo ; 1 1 ('3'p)n 1 1 oo ) ||-^^p||oo) ||Cp||oo) ||-So-||oo) ||Co-||oo "^Cy/L, 

A major difference from the existence part of the proof is the absence of 
cross-terms in the energy identities (since a; = x in the notation of Chap- 
ter [2]). In addition to that, we work in a lower order energy space and we can 
thus use the above uniform estimates to bound the term [fn]+ by Ca/L in 
L°°-norm. This observation is crucial when estimating h*. Knowing that the 
e-dependence comes only from the estimates of the cross-terms (cf. f l4.109p . 
f l4.110l) . dHUD, dmSD, fl4:T2TD and 04.1221) ). we conclude that the con- 
stants on the right-hand side of fl5.172p do not depend on e. Choosing L suit- 
ably small fl^TT^]) implies f^g* <CS* implying S*{t) <C J^S*{s)ds, since 
S*{0) = 0. By Gronwall's inequality, we conclude S*(t) = 0. In addition to 
this the conservation law (15.1670 gives ||t||2<C£^*. This estimate follows in 
the same way as (I5.166p . Thus {u,p) = {v,a). This finishes the proof of the 
uniqueness claim. 

Continuity. Integrating the identity (13.680 over the time interval [s,t], 
we obtain 

Js Js Jn Js JT"-i 

(5.173) 
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However, since {u\p^) {u,p) strongly in the energy space, we may pass to 
the hmit in fl5.173p to conclude 

S,{t)-S,{s)+ [ V,{T)dT= [ [p + R-[ [ Q + S + f. (5.174) 
Js Js Jn Js Jt"-i 

Here ^ = Zl|^|+2s<2fc-P('9^"'P'/M.s)' where is defined by dropping the in- 
dex m in the definition fl3.63p of fj^^. The terms Q, R, S and T are defined 
analogously. We claim that 




P + R- 



s Jn 




Q + S + f 



<C / ^/£,{T)V,{T)dT. (5.175) 



The inequality follows easily from the energy estimates in Chapter |H We 
observe that the estimates involving e on the right-hand side, are used only 
when estimating the cross-terms (cf. fl4.1U9l) . fl4.110p . 04.1181) . fl4.119p . 
fl4.12ip and fl4.122p ). However, the cross-terms vanish as m goes to oo (since 
X = uj = dj!p). As a result, we obtain the estimate (I5.175P with the constant C 
on the right-hand side which does not depend on e. Using ( ]5.174p and ( ]5.175p . 
we obtain 

S,{t)-S,{s)+ I V,{T)dT <C I ^fSjT)V,{T)dT. (5.176) 
J s J s 

In addition to that, for any < s < t < we have 

\£e{t)-£e{s)\<C [ V,{T)dT fl+ sup \/sJs)) — ^0 as s^t, 

Js ^ 0<s<t^ ' 

since supo<s<tE -\/^e(s) < vT. This finishes the proof of Theorem 15.11 □ 



6 Global stability 

Proof of Theorem M.^ We exploit the estimate (15.1760 to prove the theorem. 
We shall abbreviate {S,V){u,p]p){t) =:{S,V)(t) and {S,,V,){u\p^]p'){t) =: 

Existence. Let M <L/2 where L is given in Lemma [4.11 Let {u'^,p'^) be 
the associated solution to the regularized Stefan problem on the time interval 
[0,t^] given by Theorem 15.11 Set 

T: = sup|t: sup S,{u',p'){s) + [ V,{u%p'){T)dT <2m] . 

t l- 0<s<t Jo J 
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Theorem 15.11 guarantees T>t^>0. For any t<T, the estimate (15.1761) with 
s = imphes 

S,{t)+ [ V,{T)dr<S,{ul,pl) + C sup Shs) I V,{r)dT, (6.177) 
Jo o<s<r Jo 

and thus 

sup £,{t)+ [ V,{T)dT<M + cV2M [ V,{T)dT. (6.178) 

0<t<T Jo Jo 

Choose M <mm{-^^,L/2}. Inequahty (16.1781) imphes 



f 4 

sup £,{t) + / V,{t) dr < -M < 2M, 
<t<r Jo ^ 



o<i<r 

which would contradict the choice of T in case T were finite. Thus T=oo 
and the estimate (ll.30p follows easily from (16.1771) and the above choice of 
M. This proves the theorem. 

Uniqueness. We want to prove uniqueness in the class of functions 
(m,p) satisfying supo<t<oo^e('"5P)(^) + Jo^^eiu, p){r) dr <2M , where M may 
be chosen smaller if necessary. It is done in exactly the same way as the 
uniqueness proof in Theorem 15. 1[ □ 

Proof of Theorem li.it Claim 1 : Let i^' < M be any positive number, where 
M is given by Theorem II. 2[ If the initial data (mo,Po) satisfy 



S{uo,po)+ I Po- uo(l + 0'po) 



K 

then there exists a unique global solution to the Stefan problem (II. lip - 

(I1.16p . Moreover, we obtain the global bound supo<t<oo^(^) + Jo'^{T)dT <K. 
Proof of Claim 1. Let {(m'^,^'^)}^ be a family of solutions of the regular- 
ized Stefan problem satisfying the given initial condition (u''(a;,0),p''(x',0)) = 
(mq,Po)) where we choose (Mg,Po) so that 

K'Po)-^(«o,Po), ^{uI,pI]pI)->S{uo,Po]Po) as e^O 

and 



J] el \V^d^,Apl\<V~e. 



|/x|+2s<2fc 
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<f. The- 



Thus for e small, we have £e{uQ,pQ)+ J^„_iPq — J^Mq(1 + 0'pq) 

orem [L2] guarantees global existence of the solution {u'^,p'^) and also gives 
the estimate suPq^^^^ S^(t) + j^'D^{T)dT <K. Since £<£e and 'D<V^, we 
obtain supQ^-i.^^£{u^,p^)(t) + f^V{u'',p''){T)dT<K. Passing to the limit as 
e— >0, we obtain the solution {u,p) to the original Stefan problem (11.111) 
- (11.161) . The uniqueness claim follows by setting e = in the proof of the 
uniqueness statement of Theorem II. 2[ This finishes the proof of Claim 1. In 
the same way as we derived the inequality (I6.177p . we deduce for any t > 0: 

sup S,{t)+[ V,{T)dT<S,{ulpl)+CVK [ V,{T)dT. (6.179) 
0<T<t Jo Jo 



If we choose K < min{ ^^^^ ,M} =: Mi, absorb the right-most term into the 
left-hand side and drop the supremum sign, we obtain for any t>0 



We let e — > and by lower semicontinuity and the assumptions on initial data, 
we obtain 

1 



em- 



V{s)ds<S{uo,po). 







In addition to this, we obtain the conservation law 



dt{ [ p} = d,{ [ u{i + 4>'p)}. 



(6.180) 



(6.181) 



Let us set M*:=^, where Mi is defined in the line after fl6.179p . 
and assume S{uo,po)+ /]p„_i po — /q^o(1 + 0'po) <M*. Claim 1 guarantees 
the global existence of the solution (m,p) and also gives the global bound 
s\xpQ^^^^8{t) + j^V{T)dT <^. In order to prove (11.251) . we first fix any 
s > 0. The idea is to solve the Stefan problem with the new initial data 
(M^(x,0),p^(a;',0)) = (u(a;,s),p(a;',s)). The problem allows for unique solu- 
tions by Claim 1, since 



«o(l + 0Vo) 



■■£{u,p){s) + 



PO- / Mo(l + 0Po) 

Jn 



Ml Ml _ Ml 
-~^T2"~^' 
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In addition to this we have the global bound supQ^f.^^£{u^,p^)(t) + 
/p V{u^,p^){T)dT <Mi (again by Claim 1). We are thus in the unique- 
ness regime and we conclude (u^,p^)(t) = {u,p){t + s) for any t >0. We may 
now use the estimate (16.1801) to obtain (ll.25p . 

The second main ingredient in proving the decay is to control the instant 
energy in terms of the dissipation, i.e. to prove that there exists a constant 
C > such that £{t) < CV{t). We know that for +2s < 2fc 

\\d''^yp\\m<CVv. (6.182) 

Thus, the only non-trivial term left to estimate is ||'u||i2(Q). 

Claim 2: There exists a constant C>0 such that ||'u||L2(n) <C\fV. 
Proof of Claim 2. Let x G and x' G T"^^ be arbitrarily chosen. By the mean 
value theorem 

u{x)=u{x')+ / Vu{tx+{l — t)x')-{x — x')dt. 
Jo 

Note that j^n-i u{x') dx' = because m = V • (j^^ on J"--^, We thus integrate 
with respect to x' over T"~^ and then with respect to x over Q to obtain 

|T"-^| [ u{x)dx= [ [ [ Vu{tx + {l-t)x'){x-x')dtdx'dx. 
Jn JnJr^-T-Jo 

Therefore 

1 



u{x)dx 



n 



< . \, max \x-x'\\n\\T'-'^\\Wu 



|L°°(n) 



< C| I Vn| |ioo(f,) < C\ I Vu| < CVV. 

By the Poincare inequality, 

\\u\\LHn)<C\ \ / u\\L2^^)+C\\Vu\\L2(n)<CVv. 
Jn 

This finishes the proof of Claim 2. As explained above. Claim 2 and the 
estimate (I6.182p together, imply that there exists C > such that 

S<CV. (6.183) 
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Plugging 06.1831) into fll.25p yields for any s>0 and some constant a>0: 

S{t)+a I S{T)dT<S{s). (6.184) 

J s 

As in [15], p. 135, define a function V{s):= f°^S{T)dT. From fl6.184p . 
aV{s)<S{s), 

V'{s) = -£{s)<-aV{s) 

and thus V{s) <V{0)e~^°'. We integrate fl6.184p with respect to s over the 
time interval [t/2,t] to get 

Thus £{t) < ^e~t*. There exist ki,K2 > such that for any t > 

S{t)<he~^^K (6.185) 

Integrating the conservation law fl6.18ip implies Jj„_i{p{x' ,t) — p)dx' = 
J^u{l + (j)' p) . By an argument analogous to (15.1660 . \\p{t) — p\\2<CS{u,p){t). 
Combining this inequality with (I6.185p . we conclude 



£{u,p){t) + \\p{t)-p\\l<K,e 



-K2t 



for some new constant Ki > 0, K2 as in (I6.185P and for all t>0. This finishes 
the proof of the theorem. □ 
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